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Admiralty Computing Service’ 


In 1942, in order to use more efficiently the scientific staff available in the 
Admiralty, the Director of Scientific Research set up, within the branch 
directed by Dr. J. A. CARROLL, an Admiralty Computing Service to cen- 
tralise, where possible, the computational and mathematical work arising 
in Admiralty Experimental Establishments. 

Mr. JOHN Topp undertook the organisation and supervision of the Serv- 
ice. By agreement with the Astronomer Royal additional staff were attached 
to H. M. Nautical Almanac Office to carry out the computational work under 
the direction of the Superintendent, Mr. D. H. SapLer. In addition, ar- 
rangements were made to permit the employment of experts from the Uni- 
versities and elsewhere as consultants. 

The work undertaken by Admiralty Computing Service was in general 
of one of two classes: heavy computation, or difficult mathematics. Alto- 
gether more than one hundred separate investigations were carried out 
ranging from projects involving several thousand hours’ computing to small 
problems for which a solution could be obtained in a few hours. In addition, 
a considerable amount of advisory work has been undertaken, usually in- 
formally. For instance the Five-figure Logarithm Tables,reviewed in RMT 188, 
were designed by Admiralty Computing Service for the Ministry of Supply, 
as one item in a comprehensive program for providing the optical industry 
with the tables they required. For various reasons, mainly owing to the in- 
creased use being made of machines and to the availability of the U. S. re- 
print of PeTERs’ seven-figure table of natural trigonometric functions, other 


| tables were never published; but see under UMT 57. 


Shortly after the formation of Admiralty Computing Service it became 
apparent that research work in Admiralty (and other) Establishments would 
be greatly facilitated if their members were informed in certain mathematical 
and computational techniques not usually covered in undergraduate courses, 
and of which no adequate account was available in easily accessible litera- 
ture. Accordingly the preparation of a series of monographs of an expository 
nature was begun [see ACS 53, 68, 71, 101, 102, 106 (revised edition of 
ACS 26), 107]. 

It was soon realised that while centralisation within a Department was 
an improvement, nothing less than centralisation on a national scale could 
be really efficient. At the end of 1943, an approach was therefore made to 
Sir Epwarp V. APPLETON, Secretary of the Department of Scientific and 
Industrial Research, asking for consideration of the formation of a National 
Mathematical Laboratory. Discussions, in which the experience gained by 
Admiralty Computing Service played an important part, have now resulted 
in the formation of a Mathematics Division of the National Physical Labora- 
tory. Staff have been released from Admiralty Computing Service to form a 
nucleus for the computational sub-division of the new organisation. It is 
anticipated that the computational needs of the Admiralty will be met by 
outside organisations and the mathematical needs by an even larger use 
of the service of consultants, working under the general direction of the 
Director of Physical Research, Admiralty. 
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Among the consultants employed during the war were Dr. N. ARoNszajn, 
Professor W. G. BicKLey, Dr. L. J. Comrie (Scientific Computing Service, 
Ltd.), Professor E. T. Copson, Dr. J. Cossar, Dr. A. Erpéty1, Professor 
P. P. Ewatp, Dr. H. Koper, Dr. J. MARSHALL, Dr. J. C. P. MILLER, 
Professor E. H. NEVILLE. 

It is perhaps as well to explain that Admiralty Computing Service started 
on a scale severely limited by difficulties of staff recruitment; its main object 
in its early stages was to obtain the numerical results required in problems 
of war research and to make those results available to the particular estab- 
lishment concerned as early as possible. Publication of reports was then a 
secondary matter, and, in fact, it has generally been so regarded as far as 
the purely computational work is concerned. For this reason, it was not 
until the beginning of 1944 that the numbering of Admiralty Computing 
Service reports was systematised in the SRE/ACS series. All reports 
issued prior to that date, either with an NAO (=Nautical Almanac Office) 
serial number or with an SRE/MA (=Headquarters) reference number, 
were renumbered in the new series and will be referred to here by those 
numbers. 

All the reports mentioned below were issued by the Admiralty Compu- 
ting Service of the Department of Scientific Research and Experiment (Ad- 
miralty), Great Britain; generally the computational reports were prepared 
and reproduced by H. M. Nautical Almanac Office, while the mathematical 
reports were edited and reproduced at Headquarters. 

The following 21 SRE/ACS Reports have already been reviewed in 
MTAC under the heading of Recent Mathematical Tables (one under 
UMT): 


ACS RMT No. p. ACS RMT No Pp. 
19 262 13 36 65 226 12 446 
21 260 13 35 68 206 11 424 
22 267 13 39 71 206 il 424 
31 260 13 35 82 252 13 31 
39 260 13 35 4 80 
46 260 13 35 91 [UMT 41] 13 52 
47 334 16 175 93 14 70 
52 268 13 97 333 16 174 
53 206 11 424 102 288 14 76 
55 260 13 35 108 352 17 215 
62 266 13 38 


There follows a summary of 19 other items of work (ACS 7, 8, 9, 18, 20, 
26, 37, 40, [47], 51, 80, 89, 96, [97], 101, 106, 107, 109, 110, 111, 112) which 
appear to have a certain permanent value to the computer and mathe- 
matician. It is hoped to arrange for the publication in full of some of the 
reports. Suggestions as to those which appear suitable for this treatment will 
be appreciated. Details of other work, not conveniently described, apart from 
its background, will appear elsewhere. 

Copies of the reports are only available for distribution to Government 
Departments and similar agencies but arrangements have been made for 
copies of some of the reports to be deposited with the Editors where they 
may be consulted. A very limited number of photostat copies of the un- 
published tables is available for distribution or loan to institutions or indi- 
viduals with a special computational requirement. 
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7. Summation of Certain Slowly Convergent Series. Stencilled typescript, 
on one side of 4 p.; undated but issued January 1943. 20 X 32.5 cm. 


This note draws attention to a device which was apparently first applied in computa- 
tional work by P. P. Ewatp in his work on crystal-structure, Ann. d. Phys., v. 64, 1921, p. 
253-287, Gesell. d. Wissen., Gottingen, math.-phys. Kl., n. s. v. 311, no. 4, 1938, p. 55-64. 
The analytical basis of the device is the Jacopi1 Imaginary Transformation of Theta-func- 
tion Theory (see E. T. WHITTAKER & G. N. Watson, A Course in Modern Analysis, fourth 
ed., 1927, p. 475); there is a physical basis, too, which consists in replacing (taking the 
electrostatic analogy) point charges by Gaussian space-distributions. Applied to the series 


S= 2, (n + $)“e-%(™+® accuracy comparable with that obtained by summation of 50 terms 


of the original series may be obtained by taking a single term of one of the two infinite 
series into which S is transformed, and three terms of the other. 


8. Mechanical Quadratures. Stencilled typescript, 11 p. + 1 diagram; un- 
dated but issued December 1942. 20 X 32.5 cm. 


Part I is an exposition of the method of (Approximate or) Mechanical Quadratures in 
which an estimate for Ja°f(x)w(x)dx is given as ZA,f(x,) where the x, are the zeros of the 
orthogonal polynomials associated with the distribution w(x) in the interval (a, b) and the 
\, are certain constants, called the Christoffel numbers. . 

Brief discussions are given of the cases when w(x) = 1, see RMT92,132; w(x) = e~*, 
t = x*, when the polynomials are the Hermite polynomials, see RMT131,250; w(x) = e= 
when the polynomials are the Laguerre polynomials, see RMT252; and w(x) = x*e* when 
the polynomials are the Sonine polynomials. 

_ Part II applies these methods to a particular case of estimating a probability integral 
of the form t = x*. 

A third part is in preparation; this will deal with the two-dimensional case and in it an 
account will be given of some recent Russian work. 

The most promising of the methods discussed appears to be the Laguerre case and con- 
sequently the definitive table of the An, xs, referred to above, see RMT252, was prepared. 


9. Table of f(x,y) = f e-* 649 Stencilled typescript, 


2 p.; undated but issued April 1943. 20.5 X 33.2 cm. 
The function f(x, y) is tabulated for x, y = [0(.25)5; 34D], without differences, with 
a reservation that “the tabular values are unlikely to be in error by more than one unit 
in the last figure retained.” 
The table was computed from the relation 
e*f(x, y) = evi(x) 
where j is regarded as an operator such that 
(2r)! 
= xt 


where J,(x) is the Bessel function of purely imaginary argument. The BAASMTC values 
of i,(x) = x-*I,(x) were used in the computation. 


I(x) 


18. Cable Tables. Stencilled typescript, 10 p.; undated, but issued in August 
1943. 33.5 X 40.5 cm. With separate introductory text, mimeographed, 
2 p. 20.5 X 33.5 cm. 


These tables are a re-issue of earlier tables prepared by H. M. Nautical Almanac Office 
prior to the formation of Admiralty Computing Service. 

The problem is connected with the form of a heavy cable in a uniform stream, but it is 
thought that the tables are of more general interest and may have other applications. 
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Defining 


sin UdU 


the quantities tabulated are: 


f(U)dU 

and the difference (s — y)/g, for the ranges: 6 = 0(1°)90°; u = .05, .1(.1).5, .4(.2)2(.5)5(1)12, 
with 3D in y/g and s/g, but 4D in (s — y)/g for 41° < @ < 90°. 

No differences have been provided and interpolation is not always linear, though 
(s — y)/g behaves smoothly even: for small » and @ near 90°, where both y/g and s/g are 
not easily interpolable. It is not expected that the last figure will be in error by more than 
one unit, though no great effort was made to ensure end-figure accuracy. 


20. Trajectories of a Body Moving with Resistance Proportional to the Square 
of the Velocity. Stencilled typescript; 6 p. + 1 p. diagrams. 20.2 33 cm. 
Tabulations are connected with projectile trajectories when the motion is under gravity, 


with a resistance proportional to the square of the velocity. 
The functions tabulated are 


sec? ydy y- tan sec? ydy S= sec? ydy 
— Fla)’ 6 F(y) — Fla) ’ % F(y) — Fla)’ 
VF(y) — Fla)’ dT — Fle)]’ 


where = sec @ tan + In (sec + tan 0) = ydy. 


In tables I(a) to II(b), @ = 89°, and since this makes most of the values of X, Y, S, T 
negative, the quantities actually tabulated are —X, — Y, —S, —T, S. The tabulations are 
generally to 3D. No attempt at great accuracy has been made in the calculations. In the 
main tables it can be stated that 


(i) the maximum error possible is ten units in the last figure retained: this error, if it 
occurs at all, will be systematic and will therefore not enter with its full weight. 
(ii) errors of more than three units in the last place are unlikely. 


Tables have been prepared for a = — 85°(5°)85° with @ as independent variable, and 
for a = — 90°(5°)—80° with T as independent variable. 
Table I(a) : a = — 10°(5°) + 85°, @ — a from 0 to 10°. 


Table I(b) : a = 70°(5°)85°, 6 — a small, less than 0°.1. 

Table II(a): a = — 85°(5°) + 80°, 6 from —5° (or some larger value depending on a) to 

89° for a < — 56° and to 85° fora > — 45°. 

Table II(b): @ from —85° to 85°, @ from 85° to 90°. 

Table III : a = — 90°(5°) — 80°, argument T = 0(.1)1(.2)2.8. 

Table IV: an auxiliary table giving F(@) for @ = wach 1)70°; 4D], and cos? 6F(@) for 
= [60°(0°.1)90°; 4D]. 

The report contains instructions for the use of the tables and indicates the method of 

computation. 


26, 106. E. T. Copson, The Asymptotic Expansion of a Function defined by 
a Definite Integral or Contour Integral. Mimeographed, 1943, 45 p. 106. 
Second revised ed., 1946, 63 p. 20.2 X 33 cm. 


This monograph gives an account of the methods used in the asymptotic evaluation of 
integrals. It includes the method of integration by parts, Laplace’s method, Kelvin’s prin- 
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ciple of stationary phase, the use of Watson’s lemma, method of steepest descent and saddle 
point method. Applications are made to the gamma function, incomplete gamma functions, 
Bessel functions, scattering of sound waves, Airy integrals, Legendre and Hermite poly- 
nomials, and other functions. 

In the second edition the section on Airy’s integral has been revised in order to bring the. 
notation in line with that used by the British Association Tables, and an application to a 
problem in probability has been added. The Bibliography (with brief notes on the scope of 
each item) will be of assistance to those who require further information. 


37. J. M. Jackson, An Electronic Differential Analyser. Mimeographed, 
1944. 19 foolscap leaves + 4 plates. Reprinted by the Navy Department, 
Washington, D. C., Office of Research and Inventions, July, 1946. 


The processes required in a differential analyser, e.g. addition (subtraction), multiplica- 
tion, integration, and differentiation, can all be performed by simple electronic circuits. 
The chief drawback is that the lower frequencies being differentiated suffer attenuation and 
phase shift, and that the process of integration can occur only for a limited time. Both these 
troubles can be appreciably relieved by liberal use of negative feed back. 

A differential analyser was built on these lines and, within the limits of its accuracy, 
about 5%, gave good service. The main advantages, apart from the low initial cost, are the 
rapidity of operation, and the simplicity of setting up, each unit, whether adder, multiplier, 
or integrator, being simply plugged into the correct position in the chain. The output of the 
machine operated a pen and ink recorder. It is felt that a more detailed investigation of the 
possibilities of the method should be made in cooperation with experienced electronic 
engineers. 


40. Rangefinder Performance Computer. Stencilled typescript, July 1944, 3 p. 
20.4 X 33 cm. 


This device calculates the mean error and root mean square error made by an operator 
being trained in the use of a mechanism such as a range-finder. It thus reduces greatly the 
labor both of selecting trainees according to natural ability and of assessing the value of 
their training. 

The true reading, and that obtained from the operator’s use of his mechanism, are fed 
simultaneously into a differentia! gear which rotates a uniselector. The instantaneous error 
appears in the uniselector expressed as an integer up to +24. At instants determined 
either by a timing device or by the operator himself, this error is transferred to the com- 
puting mechanism, which, by uniselectors and relays, adds it to the sum of all previous errors, 
and by means of a built-in table of squares, adds its square to the sum of the squares of the 
previous errors. It also counts the number of errors inserted. These operations require up to 
4 seconds; the results, including the sign, are shown on an illuminated display panel. The 
machine can accommodate up to 100 errors, or can be set to stop automatically at any 
previous number. It can be reset to zero almost instantaneously. 

The calculation is completed either by a hand machine or by inserting the numbers 
from the display panel into a suitable network of inductance potentiometers. 


[47]. Tabulation of the Function f(x,y) = f ‘CU o(kx) cosh (ky) — 1 ]dk/sinh k. 


Photo-offset print of handwriting and machine-printed tables, February 
1945. 9 p. + 1 folding diagram. 21 X 34.5 cm. See RMT 334. 
The function f(x, y) is the solution to a two-dimensional potential problem, being that 
solution in the strip 0 < y < 1 of the differential equation 
14 & 
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which satisties the conditions = 0 on y = for = (1 + on y = 1, and 


(2y + (ZY ~0(s55). both near the origin and for large values of x. 


The function is tabulated for x = 0(.1)5; y = [0(.1)1; 4D], with second differences in 
both the x- and y-directions. There is also a diagram showing the contours in the x, y-plane 
for f(x, y) = — 1(.1) + .3, and the values of f(x, y) for x = 0(.2)5; y = 0(.2)1 obtained by 
the use of relaxation technique. 

Expansion as a series gives: 


Sle, 9) = + (2m — + fat + (On + 


and this was used, in various forms, to obtain values to 6D on the lines x = 0,x = 5, y = 0 

° and y = 1 and for one or two interior points as a check. Computation for the remaining 
interior points at interval .2 in both x and y was pérformed by applying the technique of 
relaxation to solve the simultaneous equations arising from the finite-difference equivalent 
of the partial differential equation defining the function. As far as it is known this is the 
first application of the relaxation method specifically including in the equations to be solved 
the corrections for the effect of higher order differences. Their inclusion allows the use of a 
larger interval than would otherwise be permissible. The excellence of the agreement of the 
independent calculations and the method of computation suggests that no value is in error 
by much more than one unit in the last figure retained. 


51. Table of Angular Quarter Squares. Photostat, July 1944. 11 p. 20 x 25 
cm. 


This table gives, for the range 0(3’)100°, the angular equivalent in degrees and minutes 

of the radian measure of the quarter square of the radian measure of the argument. It was 

’ prepared for the Admiralty Compass Observatory for use in analysing compass deviations. 

The preparation of the table is trivial and it was actually built up from a constant 

second difference on a sexagesimal National accounting machine, the machine being ar- 
ranged to print the final copy directly. 


80. Probability Charts for Destructive Tests. Sheets mimeographed on one 
side only, November 1945. 5 p. of introductory text, 4 folding charts. 
20.2 X 33 cm. 

Tables are given to 2D of log [ 1000 fe. ad as a function of m for the 
range of values N = 50, 100, 200, 400; M = 0(1)10, 12, 15; c = 0(1)2;” < 4N. 
Values for the limiting case N = © are also tabulated as a function of n/N. 


89. Solution of Integral Equations occurring in an Aerodynamical Problem. 
Photo-offset print of handwriting and machine printed tables, July 1945. 
17 p. 20.5 X 33 cm. 


The actual tables consist of 2 p. only, on which are tabulated: 

in Table I, K’(x), G(x), G’(x), S(x), U(x) and W(x) for x = [0(.2)10; 5D]; 

in Table II, w’(x), \’(x) and f(x) for x = [0(.2)10; 5D]. These functions are defined as 
follows: 


2v2 1 
K' k) — F, 
(x) vx +2 {( + x) Ex( ) 1(k)} 


G(x) = Fi(k), G(x) = 


Vx +2 +2 
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where Ex(k) = (1 — sin® ¢)tdg, and Fi(k) = (1 — #*sin* ¢)-tdg are the elliptic 
integrals of the second and first kind respectively and k* = x/(2 + x). 
From these S, U and W are derived by means of the equations 


S(x) = K"(x) K'(x — y)S(y)dy, = G(x) G(x — 


W(x) = — U"(x) = {G(O)S(x) — G'(x)} + f G'(x — y)S(y)dy. 
For the functions in Table II, 


1 
2)E.v(k) — Fi(k)}, 
(x) * ) i(k) )} 
v2 1 
N(x) = — 


{(2x* + 4x + 1)Ei(k) — (2x + 1)Filk)}, 


x xvx+2 


and f(x) is defined by the integral equation: f(x) = w'(x) — f r(x — y)f(y)dy. 
The problems from which the two tables arise are essentially the following: 


(g +1 — 1%) : 
h(t)dt for h(z) in 


(i) To solve the integral equation g(z) = A 
terms of the general function g(z). The solution is 


h(z) = w(s) — * S(z — t)w(t)dt 


where w(s) = 2H, and v(s) = f. g(z)(s — z)~tdz and so can be obtained by quadrature 


for any given function g(z). 
(ii) To determine the function f(z) from the equations f(z) = Pe k2(z — é)h(t)dt, where 


Simple elimination leads to the equation 
exe — = = + 2))4, 


a form identical with the first problem, the general solution of which could be used to 
compute f(z). 

The main interest in the report lies in the method of numerical solution for the integral 
equations of the second kind defining S(z) and f(z). A simple direct method is developed by 
replacing the integral by an accurate quadrature. It is also shown that the following explicit 
formulae for S(z), U(z) and f(z) may be obtained by the applicatim. of the Laplace trans- 
formation: 


1 eull-2) du du 
K,%(u) + Vu Ki(u) + 


ow du 
+ u? 
K,, I; are the Bessel functions of pure imaginary argument, primes denote derivation with 
respect to the argument, «, and the first term in f(z) arises from the fact that K,’(z) has 2 
(and only 2) zeros in the complex plane: s = — .6435 + .5012i. These expressions were 
used to provide an independent check on the accuracy of the numerical solution, and it is 
difficult to conceive a check which is more satisfactory. They do not, however, appear to 
be particularly suitable for systematic numerical computation. 
This problem has had a certain amount of publicity owing to the disparity between the 
original estimate of 300 man-years, made by the Establishment concerned, and the actual 
time of 50 hours taken by Admiralty Computing Service. 
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96. Computation of an Integral occurring in the Theory of Water Waves. 
Photo-offset of handwriting and machine printed tables, September 
1945. 6 p. 20.5 K 32 cm. 


The theory demands numerical values of 
C(u) = f° cos ou sech kh de 


where o? = gk tanh kh, for various values of the parameter h, g being a constant. It is readily 
seen that a single table of 


fv) = cos vy sech x dy 


will enable C(u) to be calculated quickly and easily for any combination of h, g and u. 
f(v) has accordingly been tabulated for v = [0(.01)5(.1)10; 4D]; no differences are 
given, but first differences are always small. 
A short auxiliary table is included of sech x with argument y where y? = x tanh x; 
sech x to 4D, with 6? for y = 0(.1)3.5. 


[97]. Tables of the Incomplete Airy Integral. Photo-offset typescript and ma- 
chine printed figures, April 1946. 5 p. of tables (photo reduced machine 
figures) and 10 p. of introductory text. 20.8 X 32.6 cm. See RMT 333. 


This report contains tables of the integral 
1 
F(x, 9) == cos (xt — 


for x = — 2.5(.1) + 4.5; y = [0(.02)1; 4D]. No A given, but the table can be interpolated 
in the x-direction using second differences; accurate interpolation in the y-direction requires 
the use of fourth differences, but the use of second differences only will give rise to a maxi- 
mum error of 2 units only. The last figure tabulated should not be in error by more than one 
unit, and all the computations have been done in units of the seventh decimal, though two 
of these were often lost in the course of the integrations. 

The title of the report is obtained by considering the function in the form: 


F(x, 9) = cos (XT 


where Y = (3y)!/8 and X = x/Y. It is thus seen that, for the purpose of extending the 
tabulation of the complete integral, the present choice of variables is not ideal. 

The report contains an account of the methods of computation used and gives the 
necessary ascending and asymptotic series. The principal method of solution was the direct 
numerical integration in the x-direction of the differential equation satisfied by the function 


=z 1 
= (x — yr*). 

A detailed account is given of the application of the National accounting machine to the 
integration of second-order differential equations of this type; the method is one which will 
have great value for use with automatic digital computing machines. 

Special methods had to be developed for use for small values of y, and generally the 
function is one of surprising computational difficulty. 


101, 107, [109, 110, 111]. H. Koser, Dictionary of Conformal Representa- 
tions, Parts I-II. Mimeographed on one side of leaves, 1945, 1946. 
36 and 48 leaves. 20.2 & 33 cm. 

Each page of this dictionary is divided into two columns; that on the left gives diagrams 


or formulae relating to z-plane while that on the right gives the corresponding diagrams or 
formulae for the w = f(z) plane. 
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In Part I the cases of the linear (f(z) = az + 5) and bilinear (f(z) = (az + b)/(cz + d)) 
transformations are discussed. For instance, explicit formulae giving the actual transforma- 
tion which carries any two non-intersecting circles into two concentric circles are given. 

In Part II, Algebraic functions and 2* for real values of a are discussed. Included in this 
part is a section on JOUKOVSKI's transformation w = az + b/z and its generalisations. 

This dictionary will be completed by the issue of three further parts: 

109 Part III, Exponential Functions and some related Functions, 
110 Part IV, Schwarz-Christoffel Transformations, 
111 Part V, Higher Transcendental Functions. 


112. ALAN BAXTER (1910-1947), The Fourier Transformer. 1947. 
This is a machine of mixed electrical and mechanical design which evaluates the Fourier 
transforms C(n) ={° f(n) cos nx dx, and s(n)= f° f(n) sin nx dx of a given function 


f(x). The integrations are carried out electrically but the selection of the wave-number n is 
mechanical. The input function is followed manually and the motion translated into a pro- 
portional A.C. voltage (50 c/s) by an inductance potentiometer. A power supply of equal 
voltage is derived from a servo-operated Variac transformer. This feeds simultaneously 22 
magslip resolvers, each of which multiplies the voltage by the sine and cosine of the partic- 
ular angle at which its rotor lies. These voltages derived from the magslip resolvers, propor- 
tional to f(x) cos mx and f(x) sin mx, are integrated by modified sub-standard K.W.H. meters. 

Each magslip rotates continuously during the transit of the input function, to produce 
the appropriate angle mx, the angles being selected by a system of gear boxes. The range of 
wave numbers 7 is from #5 to 128, a figure which can be further increased if the function is 
subdivided. One traverse of the function occupies 10 minutes, and produces simultaneously 
22 cosine integrals and 21 sine integrals whose wave numbers cover a range of 4:1. The full 
range from 7x to 128 is covered, if required, by repeated following of the input function. 

Additional dispositions of the gear boxes enable the density of wave-numbers in any 
particular range to be increased at first fourfold, and then a further tenfold, allowing for 
close investigation of any particularly interesting regions of the transform. It is hoped that 
the errors in the transform will be less than 1 per cent of its peak value. 


Joun Topp 
Univ. of London, King’s College 
D. H. SADLER 
H. M. Nautical Almanac Office 
1 Other aspects of the work of Admiralty Computing Service are described in an article 


by the present authors, in Nature, v. 157, May 4, 1946, p. 571-573; see MTAC, v. 2, p. 188. 
See also A. Erpé_y1 and Jonn Topp, Nature, v. 158, 1946, p. 690; ACS 115. 
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Seven reviews of RMT are to be found in our introductory article 
“Admiralty Computing Service,” 7, 9, 18, 20, 80, 89, 96. 


400[C].—Joser KRovAK, Achtstellige Logarithmische Tafel der Zahlen. 
Osmimistné Logaritmické Tabulky Ctsel. Prague, Geographic Institute of 
the Minister of the Interior, 1940. iv, 26 p. 14.6 X 21.1 cm. 


This little pamphlet is divided into two parts. In the first part (p. 2-14) the arguments 
are the logarithmic mantissae from 0000 to 6389 corresponding to the numbers 10 000 000 
to: 43 541 160. In the second part (p. 15-26) are given the mantissae of log N, for 
N = 4340(1)10009. In columns headed d, of each part, are the greatest and least differences 
which arise in successive lines. Examples in German and Czechish illustrate the interpola- 
tion process for getting the 8-figure logarithm of any number. 
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401[D].—JosEr KRovAk, Sechsstellige Tafeln der natiirlichen Werte der 
Funktionen Sinus und Cosinus fiir Winkel in Zentesimalteilung. Sestimistné 
Tabulky ptirozenych hodnot Funkct Sint a Cosini v setinném délent. 
Prague, Geographic Institute of the Minister of the Interior, 1943?, xii, 
52 p. + two slips (in German) listing errors in the first two editions as 
well as in the current volume. 14.9 X 20.8 cm. 


From the errata slips of this undated volume it would seem as if the edition before us 
was at least the third. And from the preface of the first edition (1941) of Kfovak’s Nattirliche 
Zahlen der Funktion Cotangens, of which the second edition (1943) was reviewed in RMT 362, 
it is clear that the first edition of the six-place table appeared in 1940. 

The table gives the natural values of sine and cosine for each centesimal minute, with 
differences, each page being devoted to 19. 

On p. 51, there are ‘‘Tafeln fiir die Berechnung zehnstelliger Logarithmen der Zahlen N 
und umgekehrt,” explanations following on p. 52. RCA 


402[D].—JoseFr KRovAK, a. Zwilfstellige Tafeln der trigonometrischen Funk- 
tionen. Dvandctimistné Tabulky trigonometrick$ch Funkct. Prague, Landes- 
vermessungsamt Béhmen und Mahren, second ed., 1944. viii, 271 p. 
+ 2 loose sheets of tables. 21.8 X 29.6 cm. 

b. Koeffizienten zur Berechnung der zweiten Interpolationsglieder. Koefi- 
cienty pro vypotet druh$ch inter polaéntch Second ed., Prague, 1944. 
vi, 26 p. 21.8 X 29.6 cm. 


a. This table, which appeared originally in 1928, is a twelve-place sexagesimal table of 
sines and cosines for each 10’ of the quadrant, and of tangents, at a similar interval, 0 to 45°. 
It is the first twelve-place table of the kind. But among recent tables Andoyer’s fifteen-place 
tables of these three functions (1915-1916), throughout the quadrant, at interval 10”, 
will be recalled. 

Each page is devoted to 10’ and columns with first and second differences (d:, d2) 
are given. 

The companion five-place table b, which appeared also with the original edition, facili- 
tates interpolations. There are also two loose-sheet tables for calculating second and third 
differences: — 2(1 — z)d2/2!, 2(1 — 2)(2 — 2)d3/3!. 

The author tells us that his values of the sine and cosine functions, insofar as they were 
correct, were rounded off from the fifteen-place Rheticus-Pitiscus table of 1613, Thesaurus 
Mathematicus sive canon Sinuum ad Radium 1.00000.00000.00000 et ad dena quaeque Scrupula 
secunda Quadrantis. . . . The values of the tangents (cotangents) were derived, with the 
aid of calculating machines, from the quotients of the values for sines and cosines. All values 
for each function were tested by differences, and by comparison with the tables of Andoyer. 


©. A. 


403(D].—J. T. Peters, Sechsstellige Werte der trigonometrischen Funktionen 
von Tausendstel zu Tausendstel des Neugrades. Ninth ed., Berlin, Wich- 
mann, 1944. iv, 512 p. 18.9 & 25.3 cm. 


The preface of this work is dated “Berlin, im August 1938” and this is the year when the 
first edition was published (iv, 512 p.). A third edition with corrections appeared in 1940 
and a seventh edition in 1943. The later “‘editions” are presumably nothing but “reprints.” 
They reflect however the great popularity of the centesimal division of the quadrant in 
German war applications. 

Pages 2-501 contain the numerical values of the trigonometric functions sin, tan, cos, 
to 6D, cot to 6S, at interval 09.001 or 10 centesimal seconds. On each page are the values 
of the functions for 10 centesimal minutes. Differences and proportional parts are given 
at the foot of each page. On pages 2-21 we also find here a 6S auxiliary table w cot w for 
each 09.01 to 29,00, 
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Apart from some worked out examples and constants, the final pages (503-512) con- 
tain 2 number of tabies for changes between the sexagesimal and centesimal units, between 
centesimal units and time, and mils, and for the improvement of the Gauss-Kriiger 
projection. 

Typographically the pages are very unattractive. 

The Czechoslovakian volumes of the six-place centesimal tables of cotangents, and of 
sines and cosines, have been already referred to in RMT 362 and 401. 


404[D].—-J. T. Peters, Siebenstellige Werte der trigonometrischen Funktionen 
von Tausendstel zu Tausendstel des Neugrades. Berlin, Verlag des Reichs- 
amts fiir Landesaufnahme, 1941. viii, 544 p. 17 K 26 cm. 


This was probably the last printed volume of tables coming from Peters’ hands. The 
preface is dated March 1941; he died in the following August. The fine appearance of this 
7-place table happily contrasts with that of the 6-place table of the previous review. As in 
that table the functions are the sin, tan, cot, cos and in the first 2° w cot w is also tabulated. 
Throughout the quadrant the interval is 07.001. Proportional parts are given at the bottom of 
each page. From 2¢ on, the first two decimals of each value are printed at the head of each 
column. The main table ends on p. 521. 

Pages 525-544 are devoted to trigonometric formulae, series, solutions of equations, 
least squares, mathematical and geodetic constants, tables indicating the relations between 
sexagesimal and centesimal degrees, centesimal degrees and radians. 

Thus endeth the twentieth major mathematical table prepared under the direction of 
Peters. See MTAC, v. 1, p. 168f. The 7-place table of the trigonometric functions for each 
0°.001 in the quadrant were sent forth by Peters in 1918, and reprinted in 1930 and 1938. 
For the American edition, 1942, see MTAC, v. 1, p: 12f. 


405[D, E, V].—IA. M. Seresriisxit, “Obtekanie Krylovykh Profilef 
Proizvol’noi Formy”’ (Flow past an aerofoil of arbitrary form), Akad. N., 
SSSR., Moscow-Leningrad, Inzhenernyt Sbornik (Engineering Review), 
v. 3, no. 1, 1946, p. 105-136. 16.4 X 25.7 cm. See RMT 194, v. 1, p. 390. 


There are the following tables: © 


T. 1, p. 111, cosh y, defined by equations: x = cosh ¥ cos 6, y = sinh y sin 0, for x = 0(.2) 
-6(.1).8(.05).9, .93(.02).97(.01).99(.005)1(.01)1.05; » = 0(.01).02(.02).22. Tables 1-9 
are to 3D. Coshy = [w + (w* — x?)#]#, where w = 3(1 + x? + ¥°). 

T.2, p. 114, ¥n(@) = [4$(1 + cos @)]", for @ = 0(5°)20°(10°)180°; m = 1(1)10(5)20(10) 
40(20)80. Same ranges of 0, m in T. 4-6. 

T. 3, p. 117, ¥n(0) sin 6, for 6 = 0(5°)20°(10°)180°; m = 1(1)10(5)20. Same ranges of 0, 
in T. 7-9. 

For T. 4 and a few others we need the definition of a conjugate; as follows: If a function 


Fis expressed in the form F = >> A, cos n6/", then its conjugate is E[F] = Ansin n0/é". 
n=0 n=0 


[E.g., taking the above ¥,(@) = [4#(1 + cos 6)]", we have y2(0) = + cos 
= 3(1 +2 cos @ + cos? 6) = +2.cos@ + 20) = }(}cos0° + 2 cos 6 + cos 26); 
= 3(}sin 0° + 2 sin + }sin 26) = }sin 6 + }sin 26.] 
T. 4, p. 122, conjugate of ¥n(@). 
T. 5, p. 123, derivative of ¥n(8). 
T. 6, p. 124, derivative of conjugate of y,,(@). 
T. 7, p. 126, conjugate of y,(@) sin @. 
T. 8, p. 127, derivative of ¥n(@) sin 0. 
T. 9, p. 128, derivative of conjugate of ¥,(@) sin 6. 

There are graphs of the various functions tabulated T. 1-9. 
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T. 10, p. 129, values of coefficients Ax,» in the development = cos with the 
recurrence formula Ax, = {[" — (k,— 1)]/(m + k)}Az-1,n for k = 0(1)24; » = [1(1) 
10(5)20(10)40(20)80; 4D]. 


On p- 130, similar values of coefficients Bz, in the development y,(@) siné 


, Baa sin 6, where Bin = 2); Ban 0, Bin Aan for 
k = 1(1)10(5)20, m = [1, 2; 4D]. 
S. A. J. 


406[E].—HarRVARD CoMPUTATION LABORATORY, Fifteen-place table of e-*. 
Bureau of Ships Computation Project. Publication no. 20, July 1945. 
iv, 6 leaves. 21.5 X 27.8. Out of print. 


This is a table of e~*, for x = [0(.05)30; 15D]. C. E. Van Orstranp, Nat. Acad. Sci., 
Memoirs, v. 14, no. 5, 1921, gave in T. VI, the value of e~*, x = [0(.1)50; 33-48D]. 


407[E].—H. S. Un er, “Special values of e**, cosh (kw) and sinh (kx) to 136 
figures,” Nat. Acad. Sci., Proc., v.33, Feb. 1947, p. 34-41. 17.3 X 25.7 cm. 


T.1. e*/™, +m = 1(1)4(2)8(4)16(8)32(16)64, 96, 192. T. I-III, Sinh (x/m), 
Cosh (x/m), for +m. Various checks for accuracy are stated. 


408[E, O].—K. G. Hacstroem, “Un probléme du calcul stochastique,” 
Férsékrings Matematiska Studier tillignade Filip Lundberg, Stockholm, 
1946, p. 104-127, table p. 127. 18.5 X 24.4 cm. 


This volume of studies in actuarial mathematics is dedicated to Dr. Lundberg, former 
general manager of the Life Insurance Co., De Farenade, and Sickness Insurance Co. EIK, 
on his 70th birthday Dec. 31, 1946. There is a “table of u(z),” arranged however as a 3D 

17 31 
an off of 

2520" T4175 “ 
— +++, u(s) = 22 + 1.333328 + 1.244425 + 1.252927 + 1.293529 + --- 


table of z for u(z) = 0(.01)2.4; 2 = u1Incoshu = 4u — =e +a 


409[F].—ALBERT GLODEN, ‘‘Factorizacién de numeros de la forma x‘ + 1,” 
Euclides, Madrid, v. 5, 1945, p. 620-621. 16.6 X 24.1 cm. 


In this note the author republishes some of his factorizations of x + 1 (MTAC, v. 2, 
p. 211). Unlike the previous table, this table lists only those factorizations which have been 
completely decomposed into primes. The variable x ranges irregularly over 132 values be- 
tween 1008 and 1500. 


D, H. L. 


410[F].—ALBERT GLODEN, Table des solutions de la congruence x‘ + 1 =0 
(modp) pour 350.000 < p < 500.000. Luxembourg, author, 11 rue Jean 
Jaurés, and Paris, Centre de Documentation Universitaire, 1946, 42 p. 
21.8 X 27.4 cm. Offset print. 


This table is an extension of four previous tables of CUNNINGHAM (1 < p < 100000), 
Hoprenot (100000 < » < 200000), Gropen (200000 < » < 300000) and DeLFELp 
(300000 < p < 350000), see MTAC, v. 2, p. 71-72, 210-211. As in the previous tables, this 
table gives only the two solutions x1, x2 of 


“t= —1(modp), p= 8k-+1, 


which are less than 4p, the other two being the negatives of these. The method of construc- 
tion of the table is as follows. First, one compiles two lists of numbers of the forms a? + 5? 
and c* + 2d. In each list each prime p of the form 8k + 1 will appear exactly once. In fact 
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those numbers which appear more than once are composite. The representations 
(1) p=i+P=2+ 
lead at once to the four solutions 

x =+d(a+b)/ac (mod p). 


To reduce these values to integers modulo p involves the solution of a single linear congruence 
for ¢ 
act = 1 (mod p). 


This procedure appears to be somewhat shorter than the one earlier described by the author 
(MTAC, v. 2, p. 71-72). 

Under the heading of applications the author gives complete factorizations of nine 
numbers of the form x* + 1 (RMT 366). It may be of interest to point out that another 
application of this table is its use to rediscover the quadratic partitions (1). In fact we may 
use the tabulated values of x:, x2 to construct / and m by: 


(mod 9), 
l= m(x2+%1) (mod p). 


If we now form the sequence 7, r2, --- by Euclid’s algorithm 
p= gl +, = gti tre, 


the first two r’s less than y> are values of a, b such that p = a? + b*. The numbers c, d may 
be found in like manner from p and m. It is worth noting that the table gives a list of primes 
of the form 8k + 1 between 350000 and 500000. 

D. H. L. 


411[F].—Irvinc KapLansky & JOHN RIORDAN, ‘“‘The probléme des mé- 
nages,”’ Scripta Mathematica, v. 12, 1946, publ. Jan. 1947, p. 113-124. 
16.7 X 24.8 cm. 


This interesting paper on a problem due to Lucas,! concludes with a small table of 
“ménages numbers” ta,z for x = 0(1)m and for m = 2(1)10. These numbers may be defined 
as follows. Let : 


(1) @1, G2, 23, ***, An 


be a permutation of 1, 2, ---, . The integer a; is said to be in a forbidden place if it is in 
either the k-th or k + ist place in (1). If (1) has precisely x elements in forbidden places it 
may be said to be of type x. The ménages number w,,z is the number of permutations of 
1,2, ---, of type x. When x = 0 this number is the same as the number of ways that n 
married couples may be seated at a round table, ladies alternating with gentlemen, in such 
a way that no one sits next to one’s spouse. This fact accounts for the name given the 
function u»,2 in general. The famous function u,»,o has been tabulated by Lucas‘ for 
nm = 4(1)20. 
D. H. L. 
1 Epovarp Lucas, Théorie des Nombres, Paris, 1891, p. 495. 


412[F].—Gino Loria, “Sulla scomposizione di un intero nella somma di 
numeri poligonali,’’ Accad. Naz. Lincei, Atti, Cl. Sci. Fis. Mat. Nat., 
Rendiconti, s. 8, v. 1, 1946, p. 7-15. 18.1 X 26.8 cm. 


This note gives three tables showing all partitions of each integer » = 100 into not more 
than (i) four squares, (ii) three triangular numbers and (iii) five pentagonal numbers.* 
The number » of such partitions is also given. Tables of this kind appear never to have been 
published before this. The number » is perhaps of more interest than the actual partitions 
themselves. It is surprising to find that the tables are quite unreliable. A complete recalcula- 
tion shows the following errata. 
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Table I 


for3+3+1,read3 +3+3+41; imsertS +1+1+1; for » = 2, read = 3 
insert4+4-+1+1; for» = 3, ready =4 
insert4 = 4, ready =5 
insert 7 +2+2+41; for» = 4, ready =5 
insert 5 +5 3; for v = 3, read = 4 
delete 5 +-4+4+ 2; for v = 3, read vy = 2 
for6+2+2,read6+6+2 

delete 8 +4+1+1; for» =5, read» =4 
insert6 +6+4+1; for = 4, read v = 5 
delete 7+ at end of first line 
for8+4+2,read8+4+4 


Table 


insert 4 + 2 + 2; for vy = 2, read vy = 3 

insert 3 + 3+ 3; for v = 1, read vy = 2 

insert 5 + 2 + 2; for v = 3, read vy = 4 

insert 4 +3 +3; forS +3 +2, readS +3 +1; for» = 2, read v = 3 
for4+4+3,read4+4+42 

insert 5 +-3 + 3; for v = 2, read v = 3 
for5+4+3,read5+4+2 

for8+3+ 3, read8+3+2 

for8+4+3, read&8+4+2 


Table III 


insert4+4+3+3+4+3;forv =5, read vy =6 
insert6+3+3+3+42;forv =7, readv=8 
insert6+4+3+2+2;forv = 5, read» = 6 

insert5 4+-5 +4-+2,and6+4+3 +3; delete5+4+43+2;fory =6,readv =7 


1 For positive integers pioneer numbers are given by the formulae T = 4n(n + 1), 
pentagonal by P = 3n(3n — 1). 


413(I, L, [L]—BAASMTC, Legendre Polynomials, (Mathematical Tables, 


Part-volume A). Cambridge, Univ. Press, 1946, 42 p. 22 X 28.4cm. 8s. 6d. 
[L, M], J. C. P. Miter, The Airy Integral, giving Tables of Solutions of 
the Differential Equation y’’ = xy. (Mathematical Tables, Part-volume B). 
Cambridge, Univ. Press, 1946, 56 p. 22 X 28.4 cm. 10s. With this Part-v, 
are also the following Auxiliary Tables on Cards. 

[I], Coefficients in the Modtfied Everett Interpolation Formula, BAASMTC, 
Auxiliary Tables, no. 1. Cambridge, Univ. Press, 1946, 1 p. 19.5 X 26 
cm. 6d. each, or 5s. per dozen. 

{I], Table for Interpolation with Reduced Derivatives, BAASMTC, Aucxil- 
tary Tables, no. Il. Cambridge, Univ. Press, 1946, 2 p. 19.5 X 26 cm. 
6d. each, or 5s. per dozen. 
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Some years ago the BAASMTC considered the problem of tables whose size did not 
justify their publication as bound volumes of the Committee’s series, but which it did seem 
desirable to publish, and decided to issue them as ‘Part-volumes,’ in paper covers, with the 
aim of ultimately combining them into bound volumes of cognate tables. Considerable de- 
lay—mainly arising out of the circumstances of the war—has occurred since the inception 
of the project, but the first two Part-volumes are now issued. Apart from the (not very) 
stiff paper covers, they are uniform, in page size, typography, arrangement, etc., with the 
main series. 

Part-volume A contains tables of the Legendre polynomials P(x) as follows: 


x = 0(.01)1, m = 2(1)6, exact or 7D, with & or 52, p. A6, A7. 

x = 1(.01)6, m = 2(1)6, exact, 7 or 8S, with & or 82,, p. A8-A17. 

x = 0(.01)1, » = 7(1)9, 7D, with 8 and 4¢, p. A18, A19. 

x = 1(.01)6, » = 7(1)9, 7 or 8S, with 5 or 82, p. A20-A29. 

x = 0(.01)1, m = 10(1)12, 7D, with 4%, and 44, p. A30, A31. 

x = 1(.01)6, m = 10(1)12, 6 to 8S, with & or 62, and 6 (for x=.5-1.5), p. A32-A41. 
x = 6(.1)11, m = 2(1)6, 7 or 8S, with & or 62, p. A42. 


The tables were designed by Dr. L. J. Comriz, and calculated under his supervision 
when he was secretary of the Committee. Their nucleus, however, dates back much further, 
to the Committee’s report for 1879, which gave exact values of P,(x) for x = 0(.01)1, up 
to m = 7. Ps(x) for the same range was taken from TaLLovist! and Hayasur (5 errors in . 
Hayashi), and checked. P9(x) for this range was calculated by use of the recurrence formula, 
while the values for m = 10(1)12 were calculated from the definitions, independently of any 
previous work, by Dr. A. J. THompson, who has edited the Part-volume. 

For values of x greater than 1 the functions were computed under Comrie’s supervision. 
For the smaller values of they were built up mechanically from the constant mth differ- 
ence, while for larger values of m values obtained by use of the recurrence formula were 
checked by differencing. A short introduction by Dr. Comrie gives further details, and there 
is a page (AS) of formulae. 

The integral 


(1) W(m) = cos — mt)dt 


was introduced by Ary? in 1838, and he gave values for m = — 4(.2) + 4, later extending 
the table‘ to |m| = 5.6. The function Ai(x) tabulated in Part-volume B is related to the 
above by 


Q) = dx) =~ cos (ae + 


where = (12/z*)?. 

It can, however, readily be shown that y = Ai(x) is a solution of the differential equa- 
tion 
(3) = xy, 
and it is more satisfactory to regard the functions Ai(x) and Bi(x) tabulated in the Part- 
volume as a convenient pair of linearly independent solutions of this differential equation. 
This equation is an approximation to amy second order linear differential equation over a 
limited range not including a singularity. For if the equation be reduced to the canonical 
form 


(4) y” + I(x)y = 0, 


- then near any ordinary point X we have 


I(x) = 1(X) + — X)I'(X) 


and inserting this in (4), we recover (3) after a linear change of independent variable. It 
was this fact that led Dr. HAROLD JEFFREYS to suggest the tabulation. The work has been 
carried out under the supervision of (and a very large proportion actually done by) Dr. 
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MILLER. Miller contributes also a scholarly, but very readable, introduction, dealing with 
the history of the function, description of the tables, their computation and checking, 
methods of interpolation, and the definitions and basic theory of the functions tabulated. 
The general solution of (3) will contain two functions, of which Ai(x) is taken as one, 
and the other chosen, Bi(x), is defined by a contour integral. Power series and asymptotic 
expansions for these functions are given. For large real positive values of x, Ai(x) tends to 
zero and Bi(x) to infinity. For negative real values of x the functions oscillate. It is therefore 
convenient to write 

Ai(x) = F(x) sin x(x), Bi(x) = F(x) cos x(x), 

Ai’(x) = G(x) sin y(x), Bi'(x) = G(x) cos ¥(x). 


The auxiliary functions F, G, x, and y are slowly varying functions of x, when x is negative 
and not too small, and are very convenient for tabulation. 
The actual tables are: 


T.I, Ai(x) and Ai’(x), x = 0(.01)2, Ai(— x) and Ai’(— x) for x = 0(.01)20, 8D, with 
& or 62, p. B18-B39; 

T. II, log Ai(x), 8D and Ai’ (x)/Ai(x), 7D, for x = 0(.1)25(1)75, with 8 or 62,, p. B40-B42; 
T. III, Zeros and turning-values of Ai(x) and Ai’(x), 8D, the first fifty of each, p. B43; 

T. IV, Bi(x) and Reduced Derivatives, x = 0(.1)2.5, Bi(— x) and Reduced Derivatives, 
for x = 0(.1)10, 8 to 10D, p. B44—B46; 

T. V, Zeros and turning-values of Bi(x) and Bi’(x), the first 20 of each, 8D, p. B44; 

T. VI, log Bi(x), 8D and Bi’(x)/Bi(x), 7D, for x = 0(.1)10, p. B47; 

T. VII, Auxiliary Functions, F(x), 7D, and x(x), in degrees, 6D, G(x), 7D, and (x), in 
degrees, 6D, for x = 0(.1)2.5, with 62, and +‘, p. B49; 

F(— x) and G(— x), 8D, x(— x) and ¥(— x), in degrees, 6D, for x = 0(.1)30(1)80, 
with &, p. B50-B56. 


Provision is everywhere made for interpolation to the full accuracy of the tables. Usually 
this is done by giving second differences (6?) or modified second differences (5?,). In some 
ranges, where the latter are inadequate, we are given y‘, whose leading term is 64/1000, 
and which contains contributions from higher differences, and formulae whereby the neces- 
sary correction can be applied with the aid of coefficients supplied on Auxiliary Table I 
(see below). With Bi(x), the tabular interval is too large for interpolation by differences to 
be convenient. Here we are given reduced derivatives, and Auxiliary Table II facilitates 
their use. 

The introduction contains a short bibliography, and a graph of the functions Ai(x) 
and Bi(x), their derivatives, and related functions. Formulae are collected on p. B17 and 
B48. 

Dr. Miller is to be thanked for his labours, and congratulated upon their result and the 
manner of its presentation. 

The modified Everett formula is, when | y{ — y§| > 1, 


+ Oh) = + Of: + + + Movs + Mint; 
when — v$| < 1, 
+ Oh) = fo + Of: + + + TH + vi) 
with @ = 1 — 6, and 

82, = & — 0.18454 + 0.0380828* — 0.008305" + 0.001951 — 
100074 = 84 — 0.278275 + 0.068558 — 0.01648" + 


Auxiliary Table I gives, on one side of a card, for @ = 0(.01)1, EZ} and E? with second 
differences, 7D, and Mj, Mj and T‘, 3D, and 4, along with the formulae. 
Defining 
= = (x)/n}, 
where h is the tabular interval, Taylor's series may be written 
f(x + Oh) = f(x) + Or + + + 4... 
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and 
hf' (x + 0h) = + + 267? + 3673 + 40974 + --~ 


(on early printings, the h of the second formula was unfortunately omitted). 

Auxiliary Table II is printed on both sides of the card. One side gives, for @ = 6(.01)1, 
#, 63, 64, exact; &, 7D; 6°, 6D; and 67, 5D with one D less for x 2 .5. The other gives 28, 
3@, 40%, exact; 564, 6D; 605, 5D; 76*, 4D; and 867, 3D with one D less for x 2 .5. 

The Auxiliary Tables can be purchased separately, and are likely to be found very useful 
and convenient aids in interpolation. It is to Dr. Miller that the lion’s share of credit for 
their inception and design must go. 


Imperial College, W. G. BICKLEY 


London, S. W. 7. 


1A. H. H. Tatrovist, Finska Vetenskaps-Societeten, Acta, v. 32, no. 6, 1904. 
o— K. Hayasut, Tafeln der Besselschen, Theta-, Kugel-, und anderer Funktionen, Berlin, 
3G. B. Arry, Camb. Phil. Soc., Trans., v. 6, 1838, p. 379-402. 
4G. B. Atry, Camb. Phil. Soc., Trans., v. 8, 1849, p. 595-599. 


F. E. J. Bore, & ANDRE Cutron, Théorie Mathématique 
du Bridge a la Portée de Tous; 134 Tableaux de Probabilités avec leurs 
Modes d’Emploi. Formules Simples. Applications. Environ 4000 Prob- 
abilités. Monographies des Probabilités, fasc. 5. Paris, Gauthier-Villars, 
1940. xx, 392 p. 16 X 25 cm. 175 francs. 


The title page description shows the impracticability of attempting to list all the tables 
in this exhaustive work. These tables fall into 3 main groups: 1, a priori; 2, bidding; 3, play. 
In the first group the a priori probabilities of all distributions of the 4 suits among the 4 
players, of the suits between the 2 partnerships, of the distributions of aces, aces and kings, 
etc., and the probabilities of voids, singletons, etc., are given. In the second group similar 
tables are given except now we know the 13 cards in the bidder’s hand. The third group of 
tables mainly covers the probabilities after the dummy has been exposed and we know 26 
cards, plus derived tables to cover the cases where part of a suit has been played and we are 
interested in the distribution of the remainder. 

The basic probabilities were calculated exactly as fractions, with the aid of a Pascal 


triangle complete up to ( %6 which the authors possess, and are given to varying numbers 


of decimal places, generally about 6 or 7, but up to 15 places in some cases. Formulae for 
making the calculations are given in the text, and the use of the Tables in evaluating a hand 
or selecting the best method of play is illustrated. 

The authors point out in the first chapter that ordinary shuffling is quite apt to give 
results that are far from the random ordering which is assumed in their later calculations. 
Statistical Laboratory Manx Huser 
University of California 
Berkeley 


415(K].—H. LasproustE & Mme Y. LaBroustE, Analyse des Graphiques 
résultant de la Superposition de Sinusoides. a. Tables Numériques précedées 
d’un Exposé dela Méthode d’ Analyse par Combinaisons Linéaires d’Ordon- 
nées. vi, 205 p. 21.9 X 31.6 cm. b. Atlas de Courbes de Sélectivité. Supplé- 
ment aux Tables Numériques, 35 plates. 21 X 30.8 cm. Paris, Presses 
Universitaires de France, 1943. 350 francs. 


The problem discussed by the authors is that of harmonic analysis. Suppose that a 
function y = (x) is a sum of a finite, though unknown, number of simple oscillations. 
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How may one find the periods, the phases, and the amplitudes of the terms? The authors’ 
approach to the problems is as follows: Suppose that 


(1) y = asin (6x + ¢) +’ sin (6’x + ¢’) +a” sin (0"x + + 
if such a representation is possible, it is certainly unique. Let us consider 2m + 1 equidistant 
points x9 — m, Xo — (m — 1), «++, Xo, «++, Xo + m symmetric with respect to the point xo, 
and let y_m, ¥-(m—1), ***, Yo, ***, Ym be the corresponding ordinates of the curve. We write 
0 = 2x/n, 0 = 2x/n’,---, = 2 cos (2xm/n), a'm = 2 cos (2rm/n’), 

If we set Y, = yz + y_y, we get 

Yy = aya sin (0x9 + >) + sin + ¢’) + 
‘so that any linear combination 

Rn = Koyo + Kivi + + 

of the quantities yo, Yi, ---, Ym can be written in the form 

Rm = pma sin + + sin + + 
where 

pm = Ko + Kiar t+ +Kmam, pm = Ko + Kiay +--+ + 


It follows that (with xo replaced by x) the phases and the periods of the terms composing 
Rn are the same as the phases and the periods of the terms of y. If the K’s are so chosen that 
all the numbers p except one—say except pm—are very small, the graph of Rn/pm gives the 
first term on the right of (1). To each period 6) = 27/m corresponds a “‘selective” combina- 
tion R». The numbers K,, are the Fourier coefficients of a function large in the neighborhood 
of the point @ and small elsewhere. The simplest combinations Rm are Sm = yom +++: + Yo 
+ +++ + ¥m-1 + Ym. The corresponding multipliers of the amplitudes are then 


on + am = sin (2m +1) /sin 


Tables I and II, p. 91-142, of a give the values, to 3D, of am and om respectively for 
m = 0(1)20, 2 = .5(.01)3(.02)5(.05)10(.1)15(.2)25(.5)50(1)100(5)200(10)500; m = 21(1)40, 
nm = 4(.02)5[then as above ]500; m = 41(1)50, m = 50(1)100(5)200(10)500. Other tables give 
values of multipliers corresponding to combinations whose basic element is not yy + y-» 
but Yu — 

In b the authors investigate in great detail products (superpositions) of simple combina- 
tions and give graphs of the ratio p/p max as functions of period n. Here p is the amplitude 
multiplier and p max is its maximum. 

The title-page states that H. Labrouste is a professor in the Faculty of Sciences of the 
University of Paris and that Mme. Y. Labrouste is an associate physicist at the Institut de 
Physique du Globe, of the University of Paris. 

A. ZYGMUND 
University of Pennsylvania 


416[K].—L. W. PoLwak, assisted by C. HEILFRON, Harmonic Analysis and 
Synthesis Schedules for Three to One Hundred Equidistant Values of 
Empiric Functions. (Department of Industry and Commerce, Meteoro- 
logical Service, Geophysical Publications v. 1.) Dublin, Stationery 
Office, 1947. xxxiii, 118 p. 24.5 X 30.2 cm. £2 s2. 

This work is designed to facilitate the work of fitting the Fourier series, 
= pot picosx + p2cos 2x + --- + qisinx + g2sin 2x + 

to a set of m equidistant values. 

The formulae devised for this purpose by F. W. BEssEL more than a century ago require 


the values of the functions cos (2xm/n), sin (24m/n), for all values of m = n. In spite of the 
age of Bessel’s formulae, only brief and inadequate tables have been available until recently. 
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The present work is divided into four parts, an “Introduction,” and three tables called 
respectively “Schedules for Harmonic Analysis and Synthesis,” the ‘“‘Register” and the 
“Index.” 

The first table, namely the Schedules, provides 5 decimal values of the two functions 
given above for values of m from 3 to 100 and for values of m = n. With respect to his choice 
of 5D the author says that the functions “are given to five decimal places in spite of the 
well-known fact that three decimal places are sufficient for most geophysical investigations. 
All the values were computed to seven decimal places and I thought it desirable to retain 
five for those who require greater accuracy, especially as the presence of the additional 
figures is no handicap to the normal user of the Schedules.” The author of this review ap- 
plauds this decision. He has had occasion to prepare a similar set of tables tom = 75 wherein 
8 decimal places are retained. The needs of the astronomer, for example, require much 
greater accuracy than the needs of the geophysicist, the meteorologist, and the economist 
who are frequent users of these methods. 

The Schedules occupy 56 pages of the work, but from the repetitive character of the 
values tabulated the tables could have been given in half the space. The present generous 
use of space was adopted, however, because after each value there is given an “‘identitica- 
tion number’’ which permits the determination of the corresponding angles by reference to 
the Index. 

The Register provides us, in order of size, the values of the angles 2xm/n, which are 
given in degrees, minutes and seconds to hundredths of a second, and in radians to 6D. 
Identification numbers are also provided which by reference to the Index permit the deter- 
mination of the corresponding values of the sine and cosine. 

The Index gives the values of the angles 2xm/n in terms of the identification numbers. 
These values are identical with those printed in the Register, but are arranged differently. 
The corresponding value of the sine is recorded to 5D for each angle together with the value 
of m and n which determines the angle. 

In addition to these principal tables the work gives 12 auxiliary tables used in connection 
with illustrative examples. A bibliography is appended which contains 17 references. This 
is quite inadequate and omits reference to such modern works as those of A. HUSSMANN 
(1938), P. TeEREBEs!I (1930), and K. Srumprr (1939). See, for example, MTAC, v. 1, p. 193, 
and v. 2, p. 32. 

The sines and cosines were computed twice, once by the author and again by C. HEtL- 
FRON, the latter using E. Girrorp, Tables of the Natural Sines (see MTAC, v. 1, p. 24f). 
Also about ten per cent. of the values selected at random were checked by F. E. Dixon 
using Chambers’ Mathematical Tables (edited by J. PRYDE). 

The book is well printed and should be very easy to use in the applications for which it 
is designed. H. T. D. 


417[K, M].—Nrets Arey, On the Distribution of Relative Errors from a 
Normal Population of Errors. Danske Videnskab. Selskab, Math.-fysiske 
Meddelelser, v. 18, no. 3, 1940, 62 p. 14.8 X 24.2 cm. 


On p. 61 is a table of r in the formula 


for f = 1(1)30(5)50(10)100, 120, o, P = [.001, .01, .02, .05, .1(.1).9; 3D]. There is also a 
table of (f + 1)# to 4D. 


418[L].—HarvarD COMPUTATION LABORATORY, Six-Place Tables of (a) 
Im(x), (b) Ym(x), (€) Sm(x) — Jm—a(x), Ym(x) — Ym-a(x), (€) 
(x), (f) dH, (x) /dx. Publication no. 21, July 1945. 
v, 11 leaves. 21.5 X 27.8 cm. Out of print. 

These tables are only for x = 24, 28, 32, 36, 40 (a), (b), m = O(1)x; (c), (d), 

m = 0(1)x + 1; and (e), (f), real and imaginary parts, m = 0(1)x. 
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419[L].—_NBSMTP, Tables of Spherical Bessel Functions, volume 1, New 
York, Columbia University Press, 1947. xxviii, 375 p. 20 x 26.4 cm. 
$7.50. 


Considering that of the eleven coordinate systems for which the wave equation is 
separable, the solutions in six involve Bessel functions, and that in four the orders are half 
odd integers, it is—as Professor PHitip M. MorsE remarks in the foreword to this volume— 
surprising that extensive previous tables of these functions do not exist (see FMR, Index, 
§17.21, p. 248, or MTAC v. 1, p. 234). The surprise is not lessened by the fact that one of 
these coordinate systems is the obviously important spherical polar, nor by the fact that 
the functions are often linked with the name of Stokes. 

The combination to which the analysis of the wave equation (and others) gives rise is 
x4J,,,4(x) rather than the function J,4;(x) itself, and it is a numerical multiple of the 
former function which is tabulated in the main table in this volume, namely (#/2x)*J4.(n44)(x). 

When we note that ‘ 

(x/2x)8J;(x) = sin x/x, 


(x/2x)\J_4(x) = cos x/x, 
and that the functions of higher order are expressible in the form 
P, cos x + Qa sin x, 


where P, and Q, are polynomials in 1/x (i.e., that the ‘asymptotic’ expansions terminate 
in this case) some reason—but still hardly adequate—can be discerned for the delay in 
tabulating this class of function. 

The main table in this volume (p. 2-323) is one of 


for m = 0(1)13, and x = 0(.01)10(.1)25, expanded for m = 13 in the range x = 10(.05)10.5. 
The standard of accuracy may be loosely described as 8S (at least) for x < 10 and 7S for 
x > 10. Actually, for smaller values of x and for the smaller positive orders, there are large 
ranges where as many as 10S are given. For values of x where the functions have begun to 
oscillate, it would be more accurate to say that enough decimals are given for the above 
standards to apply to the maxima. Functions of equal and opposite order are tabulated to- 
gether, each pair occupying 23 pages. 

Interpolation is provided for by second (or modified second) differences, where these are 
adequate, and a table (p. 370-375) of the Everett coefficients is given. But for the smaller 
values of x, and especially for the higher orders, the functions vary so rapidly that (even 
modified) second differences become inadequate, and, in fact, interpolation by differences 
becomes not feasible. Hence there are many empty columns headed &* (sometimes when the 
companion function has 63, provided). The need for interpolation in these ranges is met by 
the provision of auxiliary tables (p. 326-369) giving (x/2)*-x~J,(x), with &. (In this part 
of the table we encounter many columns headed 5‘, entirely blank!) 

These tables should prove of very great value. They will enable the solutions of many 
physical and engineering problems to be explored in considerable numerical detail, where 
previously the labour would have been quite impracticable for the physicists or engineers 
concerned. The account of the methods of computation and checking, coupled with the 
reputation for accuracy which the NBSMTP has earned, enables one to use these tables 
with complete confidence. 

In view of the undoubted value of the tables, and of the amount of useful material they 
contain, it may seem ungenerous to find fault. But, in the end, nothing but the best is good 
enough! In some respects the design and arrangement of these tables falls short of the stand- 
ards which the NBSMTP has itself set, and maintained in its earlier volumes. The number 
of empty headed columns has already been mentioned. Another point is inconsistency in 
breaking up the entries into groups by spaces. In most of the tables the now usual (and for 
many reasons optimum) groups of five, with breaks every five places from the decimal point, 
are used, but in places groups of six or more are given. The arrangement of p. 321 is not 
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happy, and there seems no obvious reason why additional horizontal spaces could not have 
been provided on this (exceptional) page. Attention to such detail would have required an 
infinitesimal proportional increase in the time and labour which has been expended on this 


volume. W. G. BICKLEY 


EpitoriaL Note: The second volume of NBSMTP, Spherical Bessel Functions is to 

contain the following: 

(a) | aay of Be algal for +2» = 29(2)43, x = 0(.01)10(.1)25, and +2» = 45(2)61, 
for x = 10(.1 

(b) Tables of A,(x) = 2’T(» + 1)J,(x)/x’, for x = ag 1)10, 2» = [1(1)41(2)61; 8-9S], and 
x = 10(.1)25, 2» = [1(2)61; 7S mostly], Also for negative values of », in regions where 
(42/x)§J,(x) does not difference well. 

(c) Table of roots of J,(x) and J,’(x) over the region covered by the main tables. 

(d) Other auxiliary tables for purposes of interpolation. 


The publication of this volume is expected by September 1947. 


420[L].—P. M. Woopwarp & Mrs. A. M. Woopwarp, with the assistance 
of Miss R. HENsMAN, H. Davies, & Miss N. GAMBLE, ‘‘Four-figure tables 
of the Airy function in the complex plane,” Phil. Mag., s. 7, v. 37, 
Apr. 1946, publ. Jan. 1947, p. 236-261. 17 X 25.3 cm. See RMT 260, 
MTAC, v. 2, p. 35. 


“Prefatory Remarks: The immediate need for tables of the Airy function in the com- 
plex plane has arisen in connection with theoretical work on the propagation of electromag- 
netic waves through the earth’s atmosphere, and it was for this particular purpose that the 
present tables were computed and produced as a report at the Telecommunications Re- 
search Establishment in February, 1945. It is not claimed that the tables are comprehensive, 
but they do provide material which, so far as is known at present, is not available elsewhere. 
The introductory matter, moreover, should be of interest independently of the tables, 
as it provides general suggestions with regard to interpolation in functions of a complex 
variable. Much of the labour associated with ordinary second and fourth difference bivariate 
interpolation may be avoided if suitable use be made of the Cauchy-Riemann equations 
and if the differences tabulated be modified accordingly.” 

The functions Ai(z) and Bi(z), z = x + ty, are independent solutions of d*w/dz* = zw, 
such that Ai(z)Bi’(z) — Ai’(z)Bi(z) = 1/x. The power series 


were used for calculation of the functions by means of the relations Ai(z) = pw: — qw2, 
Bi(z) = 34(pw: + qw2), where 1/p = (— 1/3)! 38, 1/g¢ = (— 2/3)! 31. 

There are 8 tables containing the real and imaginary parts of I-IV: Ai(z), Bi(z); 
V-VIII: Ai’(s), Bi’(z). These cover the region x = — 2.4(.2) + 2.4, — y = [0(.2)2.4; 4D] 
in the lower half of the complex plane. The upper half is covered by taking the complex con- 
jugates from the tables. Beneath each tabular entry are given $(A.* — A,?) and A,*A,?. It is 
hoped that, apart from any typographical errors which may have passed undetected, the 
errors in Tables I-IV never exceed 0.6 or 0.7 and those of Tables V-VIII never exceed 0.8 in 
units of the fourth decimal. 

Postscript: Since the preparation of these tables, there have been published (‘Annals 
of the Computation Laboratory of Harvard University,” v. 2, --- 1945] eight-figure tables 
of the modified Hankel functions of order one-third'—functions which are closely related to 
the Airy functions. The tables cover the region |x + iy| < 6, at interval 0.1 in x and y. 


Extracts from the text 
1See MTAC, v. 2, p. 176f.—EpitortaL Norte. 
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421[L, M].—J. P. Kinzer & I. G. Witson, “End plate and side wall cur- 
rents in circular cylinder cavity resonator,’”’ Bell System Technical Jn., 
v. 26, Jan. 1947, “Appendix, integration of Jo’ Jn(x)dx/Jn(x),” p. 70-79; 
tables calculated by Miss F. C. Larkxey. 15 X 22.8 cm. 


On p. 73-75 are 4D tables of F,(x) = SZ Jn(t)dt/Ji(t), and of Ga(x) = e-¥», for 
x ‘= 0(.1)9.9, m = 1(1)3. Then on p. 76-79 are 4D tables of J,(x) for x = 0(.1)9.9, m = 0(1)7, 
and of J;(x), for x = 0(.1)9.9, m = 1(1)6. 


422[(L, P].—V. I. KovALENKOV, “‘Obshchee reshenie uravneniia Gel’mgol’tsa 
s uchetom vliianiia zheleza’’ [A general solution of the Helmholtz equa- 
tion taking into account the effect of iron], Akad. N., Moscow, Ofdiel 
Tekhnicheskikh Nauk, Abtomatika i Telemekhanika, Organ komiteta 
Telemekhaniki i Abtomatiki, no. 2, 1939. 16.5 & 25.1 cm. 


On p. 9 there is a table of values of 
(1) 1! . 2(2!) + 3(3!) 


for y = 0(.01)5, 4D through 3.25, 3D or 5S thereafter. 
On p. 22 is a table of values of 
(2) 
! 272i) 


for the same range, 4D through 1.34, 3D or 4S thereafter. 


+ 


423(L, S]—Guy Lansravux, ‘Calcul des figures de diffraction des pupilles 
de révolution,” Revue d'Optique, v. 26, Jan—Feb. 1947, p. 24-45. 
15.6 X 23.9 cm. ' 


On p. 43-45 is a table of An(x) = La(x) = n! (4x)-"J n(x), for m = 1(1)30, and 
x = [0(.5)5(1)15; 6D]. 


If G(x) = e~[Li(x) + (x) +--+], 

Gi(x) = cos 6[Li(x) — (@/6)L3(x) + ---] + sin — (03/24)La(x) + ---], 

Go(x) = cos O[40L2(x) — (63/24)La(x) + ---] — sin O[Li(x) — (67/6)La(x) + ---], 

and I(x) = G,*(x) + G(x). On p. 32 are tables Gi(x), Ge(x) to 5D, and of I(x) to 4D, for 

6 =+ m, + 2x, x = 0(.5)5(1)15; there are also graphs of the six functions on p. 32-33. 
There are also tables on p. 37 of: 


G(x) = e*[Li(x) + +--+ + a(x) + ---], for k = 0, to 5D, and for 
k = 1, 2, 4, to 4D; on p. 36 are graphs of the four functions. 
On p. 40 are tables of: 


1 
G2(x) = — Li(x) + 5-21 (3r)*Ls(x) + ---,and I(x) = Gi2(x) + G2*(x), for x = 0(.5)5(1)15, 
the first two to 5D and the third to 4D. On p. 41 are graphs of the functions. 

In JAHNKE & EmpeE, Tables of Functions, 1945, p. 180-188, there are tables of L,(x), 
for nm = 0(1)8, x = [0(.02)9.98; 4-5D]. But the NBSMTP tables of L.(x), see MTAC, v. 1, 
p. 363f, are for m = 0(1)20, x = [0(.1)25; 10D]; m = 0(1)12, x = [0(.01)10; 8D]. Thus the 
Lansraux tables of L,(x) are new only for m = 21(1)30, x = [0(.5)5(1)15; 6D]. 


R. C. A. 
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424[M].—NrieEts ArRLEy, On the Theory of Stochastic Processes and their 
Application to the Theory of Cosmic Radiation. Diss. Copenhagen G.E.C. 
Gads Forlag, 1943. 240 p. 17.3 K 25 cm. 


Chapter 8, p. 222-227 is entitled, “On the numerical computation of ¥(x) = Jt e‘dx,” 
t = x*, that is, Dawson’s or Poisson’s integral, tables of which we have listed, MTAC, v. 1, 
p. 322-323, 422-423, v. 2, p. 55, 185. Arley’s only references are to JAHNKE & Empe, and 
Dawson. T. 24, p. 224-225, gives ¥(x), for x = [2(.01)10; 4S]. T. 23, p. 223, is of f(x) 


1 1-3 1-3-5 
= 2xe~Yy(x) = 1+ + (my +--+, x>1, for x = [2(.01)4(.1)10; 4-5S], 
t = x*. Of this table it is stated that “‘We estimate that the figures are correct within 1 or 2 
units in the last figure.” T. 24 was calculated from T. 23 by means of the formula y(x) 
= (2x)-le'f(x), = 
4 A 


MATHEMATICAL TABLES—ERRATA 
References have been made to Errata in RMT,412 (Loria); N 74 (FMR). 


105.—G. F. Becker & C. E. VAN OrsTRAND, Hyperbolic Functions (Smith- 
sonian Mathematical Tables), Washington, 1909—fourth reprint 1931. 


In this and previous editions the formulae 83, 84-on p. XV are incorrect. They should 

read 

83. tanh“ tan u = $gd-12u 

84. tan tanh u = $gd2u 
They are given correctly in the fifth reprint 1942. The same error occurs in A. E. KENNELLY, 
Tables of Complex Hyperbolic and Circular Functions, second ed., Cambridge, Mass., 1921, 
p. (230), and in various textbooks. 

R. O. STREET 


‘Mathematics Department, 


Royal Technical College, 
Glasgow, Scotland. 


106. W. W. DuFFIELp, Logarithms, their Nature, Computation, and Uses, 
with Logarithmic Tables of Numbers and Circular Functions to Ten Places 
of Decimals. Washington, 1897. See MTAC, v. 2, p. 161-165. 


A. The fact that Duffield prepared this table by copying from VEGA Thesaurus, 1794, 
is well known; see the references in MTAC (above), and also L. J. Comrir, Br. Astron. 
Assoc., Jn., v. 36, 1926, p. 341. 

DuFFIELD’s copying from Vega included nearly 300 end-figure errors, but his work is 
held in such low esteem that no one appears to have checked anything but the end figures. 
Recently an opportunity of examining the first six presented itself when reading the proofs 
of a new 6-figure table now in press for Messrs. W. and R. Chambers. The result of a single 
comparison was as follows. 


Page Number 


495 32067 the figures 058 3318 should be overlined 
500 33411 for 8994744, read 8894744 

526 41764 for 8920871, read 8020871 

530 42680 for 2224108, read 2244108 

580 57482 for 3318704, read 5318704 

686 89330 the overline on 9973340 should be deleted 
710 9680. in number, for 7180, read 9680 

710 96826 the overline on 9919910 should be deleted 
716 98771 for 5294508, read 6294508. 
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In Science, n.s., v. 7, 1898, p. 109-111, there is a devastating editorial on this table 
entitled ‘Logarithms on the ‘Spoils System.’” It is unsigned, but from the fearless vigour 
of its style, and the fact that Suson NEwWcoMB was the mathematics member of the edi- 
torial committee,! there can be little doubt of the authorship. A few quotations, even after a 
lapse of nearly half a century, may not be without interest. ‘Anybody who knows anything 
about the subject knows that useful tables of logarithms include from four to seven places. 
The number of problems in which a table of more than seven places would be used is ex- 
tremely small, and all extension of figures over what are actually used are a nuisance and a 
real hindrance. That the United States government should suddenly print for free distribu- 
tion several thousands of copies of this compilation must create, among those who under- 
stand, a strong suspicion of a dearth of other printable material.” Then, with biting sarcasm: 
“Their arrangement might have been worse, but only by printing the numbers in one annual 
report and their logarithms in the next.”’ Of Duffield’s claim that he did not know of Vega 
when the computations were begun: “‘This great work of Vega, which every tyro in comput- 
ing knows, was published in 1794. This is more than a hundred years ago, and it is not easy 
to understand how one could seriously think of repeating such a performance without finding 
that it had already been done.”’ Then, with his tongue in his cheek: ‘The author thinks he 
has discovered some serious mistakes in Vega, but he delicately refrains from telling what 


they are.” 
<. 


1 Simon Newcomb was mathematical editor of Science 1895-1903. Eprtor. 


B. Since L. J. C. had many years ago sent to me a copy of his 1926 review of Hender- 
son’s Bibliotheca Tabularum Mathematicarum, 1926, to which he refers above, it was a 
decided oversight on my part to omit a reference to it in RMT 319. The paragraph of the 
review which is here pertinent is the following: “‘However, he did not take Duifield at his 
face value, and so made the discovery of his dishonesty in attempting to pass a copy of 
Vega as his own computation. Peters was also aware [1922] of this fraud, and it was noticed 
independently by the present writer in 1924. Henderson’s hypothesis that Duffield was 
original up to 26000 is untenable. It is far more probable that in the copy of Vega’s Thesaurus 
which Duffield used some previous owner had entered the corrections given by Lefort up to 
26000. This supposition is supported by the fact that of the five errors in Duffield before 
26000 four were not given by Lefort; three were first pointed out by Glaisher, and in one 
case Lefort had omitted the asterisk which denoted that the error was in Vega as well as 
in Vlacq.” Before commenting on this I should point out that I did not with sufficient 
clearness indicate that of the two errata lists of Lefort the first, of 1858, was a list of errors 
in Vlacq, Arithmetica Logarithmetica, 1628, upon which Vega’s table of logarithms was 
based; and the second, of 1875, was simply a list of errors in Vega’s Thesaurus. In the first 
list, however, Lefort added a star to indicate where an error in Vlacq persisted in Vega. 

L. J. C. kindly reported to me that the “five errors” in Duffield before 26000, to 
which he referred, were in connection with the numbers 10033 (Lefort with * omitted), 
11275 (Glaisher, R.A.S. Mo. Not., v. 32, p. 258), 11699 (Glaisher, idem, v. 32, p. 258), 
22312 (Lefort), 24580 (Glaisher, idem, v. 32, p. 258 and v. 34, p. 471). It is true that Lefort 
1858 does not list an error in Vega in connection with 10033 but Lefort 1875 does list such 
an error. Thus in the five Duffield errata two were listed by Lefort. With regard to the error 
associated with 11275 Glaisher remarks the logarithm of 11275 is 4.05211, 65505, 49998, 
14..., and it is a matter of indifference whether the tenth figure of the mantissa be increased 
or not. But Glaisher lists the end figures of Duffield 65506 as an “error” and the “‘correction” 
65505. That this error was also in Peters was pointed out by L. J. C. in MTE 104. Thus 
Peters listed only four Duffield errata before 26000, not five. 
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107. MAurRICE Kraitcuik, Recherches sur la Théorie des Nombres, v. 1. Paris, 
1924. 


On p. 77-80 is a table of the factors of the two Fibonacci sequences 
1, 1, 2, 3, 5,8, 13, and 1,3, 4,7, 11, 18, 
The following two errata may be noted: 


for Us; = 79-149-2221, read 73-149-2221, 
for User = 44945570212853, read 269 - 116849 - 1429913. 


This latter error appears to be due to PouLET, since the entry is attributed to him. 

Table I, p. 131-191 gives, in effect, for each prime less than 300000 the exponent e of 
2 modulo p, that is the smallest e for which 2* — 1 is divisible by p. Actually, to save space, 
the table gives y = (p — i)/e. This table is the most extensive of its kind and has been 
used to a considerable extent in connection with tables of factors of 2" + 1 and other prob- 
lems involving the binomial congruence. Immediately after its publication this table was 
compared with a set of similar tables of CUNNINGHAM & WooDALL' extending to p < 100000 
and the resulting errata of 44 items appear in Messenger Math., v. 54, 1924, p. 184 (given 
also in Guide to Tables in the Theory of Numbers. Washington, 1941, p. 155.) 

The purposes for which this table is most frequently used require information about 
primes whose exponents are comparatively small. It therefore seemed desirable to find 
independently those primes whose exponents do not exceed 2000 in the range 100000 < p 
< 300000. A comparison of these results with Kraitchik’s table yields most of the entries 
in the errata list given below.? Since this comparison involves less than 1.5% of the entries 
in Kraitchik’s table the user of this table, who is interested in exponents beyond 2000, is 
exposed to considerable risk. 


For Read For Read 
101737 4 8 165233 4 92 
102043 2 9 4165313 96 672 

2104161 60 30 194867 7 217 
106649 8 4 216217 24 168 
107857 7 14 220243 3 213 
108497 8 16 246739 3 177 
108967 39 78 247381 1 217 

3109121 8 248 247531 5 185 
111487 6 102 250867 2 1 

2114601 2 6 254039 2 142 
119929 2 114 255071 1 2 

2121081 4 0 2267481 1 2 

2127681 8 152 272959 2 938 
141023 98 14 5284689 2 216 


TIETZE notes® that on p. 191, after 297 967, for 297 671, should be 297 971. For other mis- 
prints of primes see my Guide, p. [156]. 


D. H. L. 


1A. J. C. CunnincHaM & H. J. Woopa.t “Haupt-exponents of 2,” Quart. Jn. Math. 
v. 37, 1905, p. 122-145; v. 42, 1911, p. 241-250; v. 44, p. 41-48, 1912, p. 237-242, 1913; v. 45, 
114-125. 

r. A. E. WESTERN has already noted five of these in MTAC, v. 1, p. 429. 

3 This error caused the omission of the entry: 96135601, 881 in os writer’s table of 
composite solutions of 2" = 2 (mod ), Amer. Math. Mo., v. 43, 1936, p. 351; see MTE 102. 

‘ The discovery ¢ this error leads to the following factorization into primes 

243 + 1 = 3-83-739- 165313 -8831418697 - 13194317913029593 

5 Kraitchik had here 284687 = 13-61-359, which is therefore not a prime; the value of 
¥ was also incorrect. 

® Akad. d. Wissen., Munich, Abh., n.s. Heft 55, 1944, p. 9; see RMT 369. 
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108. NBSMTP, Tables of Sine, Cosine and Exponential Integrals, v. 1, 
1940. 


P. 59, argument column, for 1.1405, read 0.1405. 


109. NBSMTP, Tables of the Exponential Function e*, 1939. See MTAC, 
v. 1, p. 438. 


P. 168, x = 1.6742, for 5.33452 58202 12879, read 5.33452 58209 12879. 
P. 304, x = .2333, for .79181, read .79191. 


UNPUBLISHED MATHEMATICAL TABLES 


56[B].—Great Britain, Admiralty Computing Service, Tables of x'/4, 
Machine printed copy prepared by and in the possession 
of H. M. Nautical Almanac Office. Compare RMT 339, MTAC, v. 2, p. 
205. 


Several requirements arose for quarter powers during the course of the computational 
work undertaken by Admiralty Computing Service at H. M. Nautical Almanac Office 
during 1943-1945. In the same period the Office was faced with the training of new staff 
with no previous computing experience. It was accordingly decided to make systematic 
tables of the four powers — 3, — 4) + 4) + 2 for a comprehensive range of argument; by 
this means considerable individual calculation for special investigations was avoided and the 
new staff provided with excellent material for elementary training in computing and 
tabulation. 

Copy has been prepared for two tables in both of which the four functions are arranged 
side by sidein the order +}, —}, +2, —? for the range x = 1(.01)10(.1)100(1)1000(10)10000. 

Table A. An accurate table with at least 7S with manuscript first differences written in 
small figures interlinearly. The number of decimals retained is: 


Range Power 
x +t +i 
1-10 6 7 6 7 
10-100 6 7 6 8 
100-1000 6 7 5 9 
1000-10000 6 7 9 


Table B. A “working”’ table to 11 or 12S intended solely to give the tabulated values to 
the greatest accuracy to which they are available; therefore no differences are provided and 
the end-figure may be in error by several units. The number of decimals (D) and the error 
(E) in the last figure which is unlikely to be exceeded are given in the following table: 


Range Power 
* +t +3 
DE DE DE DE 
1-10 10 2 10 2 10 2 10 2 
10-100 10 3 10 2 7 3 11 4 
100-1000 9 2 10 2 8 2 12 5 
1000-10000 10 2 2 i2 2 


The original aim was to provide a table giving 7S accuracy throughout, interpolable 
with only trivial second difference corrections. Basic values were calculated to 10 or 11S 
for x = 1(.01)3(.05)6.5(.1)10 for powers + 4; x = 1(.01)4(.05)7.5(.1)10 for powers + j. 
These were multiplied by the appropriate powers of 10 to give powers of 10x, 100x and 
1000x over the same ranges of x. Values of all 16 functions were then obtained for a 
uniform interval of .01 in x, over the whole range x = 1 to 10, by standard methods of 
interpolation to fifths and tenths on the National machines. The copy in each case was 
prepared by integrating on the National machine from differences produced by end-figure 
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differencing of the interpolated values. Each table contains 72 foolscap pages, each with 
4 X 50 entries, apart from differences. 
Joun Topp & D. H. SapLer 


57[D].—Great Britain, Admiralty Computing Service, Six-figure logarith- 
mic-trigonometrical tables. Copy prepared by machine printing with ms. 
proportional parts, and in the possession of H. M. Nautical Almanac 
Office. 


At the request of the Ministry of Supply, Admiralty Computing Service undertook to 
prepare copy for a six-figure table of logarithmic-trigonometric functions for use in the 
optical industry. The argument was to be in degrees and decimals. The obvious course of 
reproducing by photo-lithography Peters, Sechsstellige Logarithmen der trigonometrischen 
Funktionen von 0° bis 90° fiir jedes Tausendstel des Grades, Verlag der Preussischen Landes- 
aufnahme, Berlin 1921, was ruled out as being too expensive; it was decided to compromise 
by photographing the first five degrees of Peters’ table and adding to it a new table at 
interval 0°.01. 

It is this latter table, giving logarithms of the four functions sine, tangent, cotangent 
and cosine in the range 5°(0°.01)45°, which has been prepared. For various reasons, it was 
later decided not to proceed with the publication of the table. 

The copy was prepared by building up from the second differences of the known end- 
figures; it was prepared on the National machine in a form as closely similar to Peters’ 
table as the limitations of the machine would permit. Proportional parts were written in 
by hand. The copy of 80 p. is in perfect condition for the printer, but not good enough for 
photo-lithography. 

‘Joun Topp & D. H. SaDLer 


58[D, E].—Great Britain, Admiralty Computing Service, Tables of 
cosec? x — x~? and x~* — cosech* x. On National differencing sheets and 
in MSS. prepared by and in the possession of H. M. Nautical Almanac 
Office. 


In connection with a special investigation tables were required of the two functions 
cosec? x — x? and x? — cosech? x 


to 14D for the range x = 0(.001)1.6. They have been calculated to 17D for x = 0(.01)1.69 
and now await interpolation to tenths. 
The fundamental values have been deduced from the tables of C. E. VAN OrsTRAND, 
Nat. Acad. Sci., Memoirs, v. 14, no. 5, 1921. 
Joun Topp & D. H. SADLER 


59[L].—GreatT Britain, Admiralty Computing Service, Tables of C(m, u; x). 
Machine printed copy prepared by, and in the possession of H. M. Nau- 
tical Almanac Office. 

These tables of the hypergeometric function 
C(m, u;x) = F(i,m — p; m, x) 


have been prepared at the suggestion of Dr. W. G. BickLEy and Dr. J. C. P. MILLER, witha 
view to their use in the summation of certain slowly convergent series. Consider a series 
un, in which the ratio of consecutive terms can be written, to a good approximation for 
values of m greater than N, in the form: 


Unsi/Un = x(1 — wn + 
where yu and » are constants. 
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Then, to the same approximation, 


Unyi/Un = x(m — p)/m, where m=n+p», 
and so : 


> uny = unC(m, x). 
A good approximation to the tail of the series is thus obtained. 

The function is tabulated to 4D for 4m = 40(1)44, 80(1)84; 4u = 4(1)20; x = .8(.02) 
-9(.01)1. 

No differences are given, since it is assumed that the table will generally be used without 
interpolation. One page is devoted to each value of x and contains the two double-entry 
tables corresponding to the two series of values of m. 

The method of calculation should result in the last figures not being in error by more 
than .6, but the system of checking adopted was not.capable, throughout the whole range, 
of guaranteeing an accuracy of more than one unit. 

The general method of computation was the repeated application of the recurrence 
formulae 


C(m + 1,4; x) = oz {C(m, — 1} 
C(m, w; x) = — 1)*{(m — (m — w)(1 — x)C(m, — 1; x)} 


due regard being paid to the loss of figures inherent in this method. Initial values were ob- 
tained from the formulae 


and 


C(l,u;x) = (1 ¥1, 
and 
C(2,1;x) = —x (1 — x) 
with 
C(m, w; 1) = (m — 1)/(u — 1). 
Joun Topp & D. H. SADLER 
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In the introductory article, ‘‘Admiralty Computing Service,” there are 
reviews of publications 37, Electronic Differential Analyser; 40, Rangefinder 
Performance Computer; 112, The Fourier Tranformer. 

There is an interesting biographical sketch, and a portrait of HOwARD 
HaTHAway AIKEN (1900- ), in Current Biography, v. 8, no. 3, March 
1947. 

In Wisconsin Engineer, v. 51, Dec. 1946, p. 10-12 is an article by WALTER 
GrauaM, “Do you know your slide rule?”’, in which the author explains the 
slide rule solution of equations of the type x* = k, a* = x, tanx = kx, and 
sin x = kx. See MTAC, v. 1, p. 203, Q 8 and v. 2, p. 194, 25. 

Dr. LotHAaR ScHRUTKA, professor of mathematics in the Technische 
Hochschule, Vienna, is the author of a third much revised edition of his 
Theorie und Praxis des logarithmischen Rechenschiebers, Vienna, Deuticke, 
1943. xii, 101 p. The first edition appeared in 1911, and the second in 1929. 
There are 32 titles in the bibliography, p. 95-96, and there is a full index, 
p. 97-101. 

H. H. ArkENn & GRACE M. Hopper, ‘“The Automatic Sequence Controlled 
Calculator,” Electrical Engineering, v. 65, Aug.—Nov. 1946, p. 384-391, 
449-454, 522-528. See MTAC, v. 2, p. 185f. 
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29[Z|.—D. R. HARrTREE, Calculating Machines, Recent and Prospective 
Developments and their Impact on Mathematical Physics. Cambridge, 
University Press, 1947. 40 p. + a double plate. 12.4 X 18.1 cm. 2 shillings. 


This is a lecture given on the occasion of the author’s inauguration as Plummer Pro- 
fessor of Mathematical Physics in the University of Cambridge. Its general purpose was to 
acquaint his audience with the war time expansion of large scale computing units in the 
United States. The lecture is divided into 9 parts by the following subheadings: (1) Intro- 
duction, (2) Two classes of calculating equipment, (3) Functions of components of a digital 
machine, (4) The ENIAC, (5) The master programmer, (6) Example of the application of 
the ENIAC, (7) Prospective developments, (8) The impact of these developments on mathe- 
matical physics, (9) Conclusion. 

Most of the lecture is devoted to a discussion of the ENIAC, see MTAC, v. 2, p. 97-110. 
Part 6 is a short description of an interesting boundary layer problem which the author put 
on the ENIAC in June 1946. It consists in solving the non-linear system of three differential 
equations 


Br" =fr+ (Wy, 
with the two-point boundary conditions, 
f=k=r =0Oatx=0, kh=2,7 =Oatz = 


Part 7 isa very brief description of the Electronic Discrete Variable Calculator (EDVAC) 
type of machine. 

Part 8 is an interesting discussion of the way that the possibility of high speed large 
scale computing alters the outlook of the mathematical physicist. An ordinary system of 
linear equations becomes a problem in minimizing a quadratic form. A second order partial 
differential equation with certain boundary conditions becomes an integral equation with 
“built in” boundary conditions. To quote the author: “The facilities offered by these new 
calculating machines will at least make the formulating of physical problems in terms of 
integral equations and variation equations more familiar and may in time wean us from our 
present tendency to regard a differential equation as the basic way of formulating the 
mathematics of physical problems.” 


D. H. L. 


30[Z].—Francis J. Murray, The Theory of Mathematical Machines. New 
York, Columbia Univ. Press, 1947. vii, 116 p. 21.5 X 28 cm. Litho- 
graphed ; plastic binding $3.00. 


This work appears to be a set of lectures on certain mathematical aspects of computing 
devices. There is a wide variety of both devices and aspects. Each of its four parts is divided 
into four or five short chapters of a few pages each. There are about 200 original line draw- 
ings that add much to the interest of the volume. Unfortunately these are not numbered 
so that sometimes the reader is in doubt as to which drawing is being referred to in the text. 
The book is lithographed from typewriting and is very neat. 

The reviewer has found it difficult to give a short account of the actual contents of the 
book. There is a great deal of detail (not indicated in the chapter headings) on some topics 
and very little on others. On the whole the book is devoted to continuous or analogue de- 
vices almost entirely; there is only the briefest mention of high speed digital computers. 
There is a great deal of space devoted to electronics but only one example of an electronic 
counter, a soft tube prewar type. There is much material of an electro-mechanical nature but 
no mention is made of the possibility of using relays for computing. The mathematical dis- 
cussion ranges in depth from the identity 


= (x — — 
to Hilbert space. The four parts may be described as follows: 
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318 NOTES 


Part I, Digital machines, is disappointingly brief (12 p.) and is divided into counters, 
adders, multipliers, and ‘‘the punch card machine.” The discussion is largely devoted to 
mechanical parts used in desk calculators such as the Leibniz wheel and Napier’s bones. 
There is a description of the Hollerith card sorter. 


Part II, Continuous operators, sounds the keynote of the text: Numerical quantities 
can be represented by physical magnitudes. The magnitudes discussed range from linear 
displacements to the phase angle of alternating currents. There is a good treatment of linear 
networks. Among multiplying devices there is a description of ‘‘square”’ gears, variably 
wound potentiometers, and rectifiers. Integrators and differentiators, both mechanical and 
electrical, are treated in great detail. The rest of Part II is devoted to the theory of ampli- 
fiers, servomechanisms, selsyn units and other electrical devices and their uses in mathe- 
matical machines. 


Part III, The solution of problems, is devoted to composite machines for solving systems 
of linear equations, ordinary and partial differential equations. These machines include 
the network analyzer, differential analyzer and two electronic computers for linear equations. 
One of the latter has been designed by the author and is fully described. The mathematical 
treatment here is particularly interesting. The reader will find this material under ‘“Ad- 
justers” (p. 84-94, unfortunately the book has no index). 


Part IV, Mathematical Instruments, is concerned with planimeters, integrometers, 
harmonic analyzers and cinema-integraphs. There is a page and a half of bibliography 
arranged topically. This does not include a large number of references inserted in the text. 

The reader, whether he be interested in mathematical machines from a technical or a 
purely mathematical point of view, will find something interesting on every page. It is 
to be hoped that a second volume dealing with the theories of the many other interesting 
devices developed during the war may be eventually forthcoming. "1 


NOTES 


73. THE CHECKING OF FUNCTIONS TABULATED AT CERTAIN FRACTIONAL 
Points.—Many functions involving a parameter v, in particular Bessel 
functions J,(x), Y,(x), I,(x), K,(x), etc., besides being tabulated for + in- 
tegral values of v, as well as for y = 0, are often’given for non-integral values 
of between —1 and 1, especially for = + ?, + 3, + 3, +} and +}. 
When it is desired to perform the equivalent of a differencing check upon 
these or related functions (e.g. the zeros of these functions) considered as a 
function of » for fixed x, due to the irregular interval in », it is necessary 
to take the divided differences. For any fixed set of m »’s, it is possible to 
obtain coefficients of f, for the last, ie. (n — 1)th, divided difference which 
should vanish if the function behaves as a polynomial of the (m — 2)th 
degree in v. Thus an error ¢ in any entry f, (this includes rounding errors) 
will usually show up by being multiplied by the coefficient of f,. 

The coefficients which are given below are for three important cases 
likely to arise in practice, especially with Bessel functions: - 


(a) 7th divided difference for f, involving the 8 points » = + 3, + 3, 
+ 3, +}. 
(b) 10th divided difference for f, involving the 11 points y = +1, + 3, 


+ 3, +4, +}, 0. 
(c) 10th divided difference for f, involving the 11 points vy = + 3, + 3, 
+ 4, + 3, +}, 0. 
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In case (c), omitting the f+3 and f+3 leaves 7 points at the uniform interval 
of } in », which might be amenable to an ordinary differencing check. It 
is too cumbersome to work with more than 11 points f,. At any rate, if to 
(c) there were added f, for »y = + 1, the set of v’s would include + 1, + 3, 
+ 4, + } and 0, which are 9 points, again at the uniform interval of }, and 
to which an ordinary differencing check could be applied. The joint use of 
any two of (a), (b), or (c), when possible, lessens the likelihood of passing 
a double error which, by compensation, might yield a small divided differ- 
ence in one case. 


Divided Difference Formulae: 
(a) » = + 2, + 3, + 3, + 3; 7th divided difference 
= [21 28896 (f: — f—-2) — 42 45696 (fi — 
+174 49344 (fy — f-i) — 183 30624 (f; — f-3)]. 
(b) » = +1, + 3, + 3, + 3, + }, 0; 10th divided difference 
T5495 [10 50192 (fs + fr) — 324 40320 (f-2 + fr) + 573 16896 


(f-1 + fi) — 2944 57680 + fi) 
+ 3910 53312 (f-: + fz) — 2450 44800 fo]. 


(c) »= +i, + §, + 3, + 3, + 3, 0; 10th divided difference 
1 
= [454 16448 (f-2 + fa) — 1637 62560 (f-a + fi) 


+ 5376 98304 (f-1 + fi) — 18845 29152 (f-a + fi) 
+ 19552 66560 (f-: +f) — 9801 79200 fo]. 
HERBERT E. SALZER 


NBSMTP 


74. Poisson’s OR Dawson’s INTEGRAL AND ANOTHER INTEGRAL.—In 
W. O. ScoumANN, Elektrische Durchbruchfeldstiérke von Gasen, Berlin, 1923, 
p. 243 and 241, are tables of f(x) = /(7 e*dt, 2-5S, z = #, of log f(x), 45S, 
and of e-“f(x), 3—6D, w = x*, each for x = [.1(.1)2.6(.2)7.4]. The refer- 
ence to Schumann in FMR Index, p. 219, cannot be verified. See MTAC, 
v. 2, p. 55, N45; and RMT 378 and 424, ArLEy. On p. 242 are graphs of f(x), 
and on p. 234-235 is a small 4D table of MS(x) = fi e~“dt, u = t*, for 
x = 0(.1)1, © ; graphsof S(x), and of MS(x),aregivenon p. 237. M = I'(1.25) 
is the value of the integral when x = ©, and is approximately .9064. More 
generally BIERENS DE Haan, Nouvelles Tables d’ Intégrales Définies. Leyden, 
1867, and New York, 1939, Table 26(4), gives {0*% edt = n“I'(n-), if 
v=". But = x") = xM(n, — x") where 
I'(p + 1, x) = e-*x7dx, or as in K. Pearson, Tables of the Incomplete 
I-function, London, 1922, y, = n[I'(1/n)}° edt = I(u, p) = IL x", 
(1 — n)/n]. For n = 4, there is a table, p. 118-126, of ys for u = 2x‘ 
= [.1(.1)27; 7D], p = — .75; also, p. 164, wu = [0(.1)6; 5D]. See also the 
little tables for this particular case by F. Emp, Z. f. angew. Math. u. Mech., 
v. 14, 1934, p. 336-339, S(x), x = [.8, .9(.01)1, 1.1; 6D]. 


Murtan S. Corrincton & R. C. A. 
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CORRIGENDA 


QUERIES 


22. INTEGRAL EVALUATIONS.—What methods are available for evaluat- 
ing the integrals 


cos (0s + axe + ass + -++)dx, and 
0 


f cos (@1x + do + + a_ox? + ---)dx/x", 
t 


except for point by point numerical integration? 


Murtan S. CORRINGTON 
Radio Corporation of America 
Camden, N. J. 


QUERIES—REPLIES 


29. TABLEs oF N?/? (Q5, v. 1, p. 131; QR8, p. 204, 11, p. 336, 13, p. 375, 
14, p. 407).—In L. Potin, Formules et Tables Numériques, Paris, 1925, there 
is a table (p. 416-417) of N*”*, N = cos 0, for @ = [0(30’)90°; 4D]. 


CORRIGENDA 


As the result of recomputation Dr. J. W. WRENCH, Jr. requests that the last three 
decimal places of each of four values given MTAC, v. 1, p. 298, 1. 11-12, be corrected to 
read as follows: ber 15,— 535; bei 15,— 887; ber’ 15,— 317; bei’ 15,— 368. 

V. I, p. 33, 472, for 1778, read 1777. P. 468, for Everett, read Everett, J. D. 

V. 2, p. 65, for Block, read Bloch; for Brendle, read Brendel. 

P. 250, RMT 362, col. “Diff. for 1,” the second, fourth, and fifth entries should respec- 
tively read: 63 to 2.9, 0.3 to 1.6, 1.6 to 15.8. These corrections are due to slips made in 


the editorial office; the ms. of L. J. C. was faultless. P. 271, 1. — 3, end of line, for 
(3), 1; read (4)1;. P. 277, 1. 6, for 54, read 55; 1. 7, for 55, read 54. 


Mr. D. F. FerGuson has now (May 23, 1947) carried on his calculation of the value of 
® to 750D, and discovered errors in Dr. Wrench’s computations. The value of # on p. 245 
has to be amended 723-743D, and the value of tan~ } on p. 247, 725-743D. No announce- 
ment will be made in MTAC of the exact corrections here necessary until Mr. Ferguson has 
completely checked the remaining 58 decimal places of our published value of x. 
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